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Abstract 

This paper extends the proposal of [Phys. Rev. Lett. 101, 203901-4 
(2008)] to invisibility carpets for rigid planes and cylinders in the context of 
pressure acoustic waves propagating in a compressible fluid. Carpets under 
consideration here do not touch the ground: they levitate in mid-air (or float 
in mid- water), which leads to approximate cloaking for an object hidden un- 
derneath, or touching either sides of a square cylinder on, or over, the ground. 
The tentlike carpets attached to the sides of a square cylinder illustrate how 
the notion of a carpet on a wall naturally generalizes to sides of other small 
compact objects. We then extend the concept of flying carpets to circular 
cylinders and show that one can hide any type of defects under such circular 
carpets, and yet they still scatter waves just like a smaller cylinder on its 
own. Interestingly, all these carpets are described by non-singular acoustic 
parameters. To exemplify this important aspect, we propose a multi-layered 
carpet consisting of isotropic homogeneous fluids with constant bulk modulus 
and varying density which works over a finite range of wavelengths. 

Pacs: 00.00, 20.00, 42.10, 42.79.-e; 02.40.-k; 41.20.-q 

Keywords: (000.3860) Mathematical methods in physics; (260.2110) Electromagnetic the- 
ory; (160.3918) Metamaterials; (160.1190) Anisotropic optical materials; invisibility; cloak. 
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1 Introduction 



There is currently a keen interest in electromagnetic metamaterials within 
which very unusual phenomena such as negative refraction and focussing 
effects involving the near field can occur [H El El H]. A circular cylinder 
coated with a negative refractive index displays anomalous resonances [5] 
and can even cloak a set of dipoles located in its close neighborhood [6] . The 
dielectric cylinder itself can be made transparent with a plasmonic coating 

m. 

However, cloaking of arbitrarily sized objects requires anisotropic hetero- 
geneous media designed using transformation optics [HI IB]- The first exper- 
imental demonstration of an invisibility cloak was obtained at 8.5GHz and 
fueled the interest in this new field of optics [10]. Electromagnetic cloaks 
also allow for mirage effects [11] and their efficiency very much depend upon 
the smoothness of their boundaries |12| . Importantly, mathematicians have 
also proposed some models of cloaks, but in the context of inverse problems 
in tomography [131 [H] and also further explained what are the appropriate 
boundary conditions on the inner boundary of cloaks [151 EE]- These latter 
works open new vistas in acoustic cloaks, and thus attract growing attention 
in the physics and mathematics communities. Whereas cloaking of pressure 
waves in two-dimensional [TBI E] and three-dimensional [201 |2T] fluids, anti- 
plane shear waves in cylindrical bodies [22], and flexural waves in thin-elastic 
plates [231 121] is well understood by now, cloaking of in-plane coupled shear 
and pressure elastic waves still remains elusive [25] as the Navier equa- 
tions do not retain their form under geometric transforms. Such cloaks might 
involve (complex) pentamode materials such as proposed in [27]. However, 
it has been know for over ten years that coated cylinders might become neu- 
tral in the elastostatic limit [2B], and this route might well be worthwhile 
pursuing more actively with e.g. pre-stressed coated elastic cylinders, as this 
would allow for a simple experimental setup. 

During the same decade, some theoretical and experimental progress has 
been made towards a better understanding of band spectra for linear surface 
water waves propagating in arrays of rigid cylinders [BU [2H1 EDI Ell E21 E3] , 
or over a bottom with periodically drilled holes. Focussing effect of surface 
water waves was also investigated by a handful of research groups in arrays of 
circular and square cylindrical holes [331 ES] , as well as using metamaterials 
designed as fluid networks [BE], [B7J- However, a further control of surface 
water waves can be obtained via an alternative route, that of transforma- 
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tion acoustics. It has been actually demonstrated that broadband cloaking 
of surface water waves can be achieved with a structured cloak, with an 
experimental confirmation at 10 Hz[38j. 

In this paper, we focus our analysis on cloaking of pressure waves with 
carpets [39J which have recently led to experiments in the electromagnetic 
context [Ml Here, we would like to render e.g. pipelines lying at the 
bottom of the sea or floating in mid-water undetectable for a boat sonar. 
These pipelines are considered to be infinitely long straight cylinders with a 
cross-section which is of circular or square shape. A pressure wave incident 
from above (the surface of the sea) hits the pipeline, so that the reflected 
wave reveals its presence to the sonar boat. We would like to show that we 
can hide the pipeline under a cylindrical carpet (a metafluid) so that the 
sonar only detects the wave reflected by the bottom of the sea. 

2 Governing equations for pressure and trans- 
verse electric waves 

For an inviscid fluid with zero shear modulus, the linearized equations of 
state for small amplitude perturbations from conservation of momentum, 
conservation of mass, and linear relationship between pressure and density 
are 



where p is the scalar pressure, v is the vector fluid velocity, p is the unper- 
turbed fluid mass density (a mass in kilograms per unit volume in meters 
cube), and A is the fluid bulk modulus (i.e. it measures the substance's resis- 
tance to uniform compression and is defined as the pressure increase needed 
to cause a given relative decrease in volume, with physical unit in Pascal). 
This set of equations admits the usual compressional wave solutions in which 
fluid motion is parallel to the wavevector. 

In cylindrical coordinates with z invariance, and letting the mass density 
be anisotropic but diagonal in these coordinates, the time harmonic acoustic 
equations of state simplify to (the exp(— jut) convention is used throughout) 
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where u is the angular pressure wave frequency (measured in radians per 
unit second). Importantly, this equation is supplied with Neumann boundary 
conditions on the boundary of rigid defects (no flow condition). 

Cummer and Schurig have shown that this equation holds for anisotropic 
heterogeneous fluids in cylindrical geometries [18] , and they have drawn some 
comparisons with transverse electromagnetic waves. In the transverse electric 
polarization (longitudinal magnetic field parallel to the cylinder's axis): 

V ■ (e; l VH z ) + cu 2 e ol i H z = . (3) 

where H z is the longitudinal (only non-zero) component of the magnetic field, 
e r is the dielectric relative permittivity, e p is the inverse of the square 
velocity of light in vacuum, and us is the angular transverse electric wave 
frequency (measured in radians per unit second). Importantly, this equation 
is supplied with Neumann boundary conditions on the boundary of infinite 
conducting defects. 

In this paper, we look at such 'acoustic' models, for the case of flying 
carpets which are associated with geometric transforms. While we report 
computations for a pressure field, results apply mutatis mutandis to trans- 
verse electric waves making the changes of variables 

p\ — ► H z , p i — ► e r , A -1 i — ► e fi , (4) 

in ([2]). However, after geometric transform, (j2J) will involve an anisotropic 
(heterogeneous) density p and a varying (scalar) bulk modulus A, see ffT2j) . 

whereas flSJ) would involve an anisotropic (heterogeneous) permittivity e and 
a varying (scalar) permeability /i, see ( II Op . It is nevertheless possible to work 
with a reduced set of parameters, to avoid a varying A in acoustics (resp. ji 
in optics), as we shall see in the last section of the paper. 

3 From transformation optics to transforma- 
tion acoustics 

Let us consider a map from a co-ordinate system {u, v, w} to the co-ordinate 
system {x,y,z} given by the transformation characterized by x(u,v,w), 
y(u, v, w) and z(u, v, w). 

This change of co-ordinates is characterized by the transformation of the 
differentials through the Jacobian: 
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dx \ 


/ du 


dy 


1 = J xu dv 


dz 


\ dw 



, with J xu = kt J - J —s • (5) 
o{u, v, w) 

In electromagnetics, this change of coordinates amounts to replacing the 
different materials (often homogeneous and isotropic, which corresponds to 
the case of scalar piecewise constant permittivity and permeability) by equiv- 
alent inhomogeneous anisotropic materials described by a transformation 
matrix T (metric tensor). The idea underpinning acoustic invisibility [9] 
is that newly discovered metamaterials should enable control of the pressure 
waves by mimicking the heterogeneous anisotropic nature of T with e.g. an 
anisotropic density, in a way similar to what was recently achieved with the 
permeability and permeability tensors in the microwave regime in the context 
of electromagnetism [TO] . 

In the sequel, we propose an alternative derivation of the analogies be- 
tween transformation optics and acoustics first drawn by Cummer and Schurig 
[T8] . Our proof involves a lemma first established in [19] in the context of 
duality relations for the Maxwell system in checkerboards. On a geometric 
point of view, the matrix T = J T J/ det(J) is a representation of the metric 
tensor. The only thing to do in the transformed coordinates is to replace the 
materials (dielectric, homogeneous and isotropic) by equivalent ones which 
now exhibit some magnetism, which are heterogeneous (dependance upon 
u,v,w co-ordinates) and anisotropic. Their properties are given by [UJ 

£ = e.T- 1 , and // = T" 1 . (6) 

In transverse electric polarisation, the Maxwell operator in the transformed 
coordinates writes as 

V x (gfV x H,) - fioeoco 2 ^ = (7) 

where H/ = H z (x,y)e z , e[ and // are defined by (jSJ). 

We would like to deduce the expression of ([3]) in the transformed coordi- 
nates from the vector equation (J7|). For this, we need the following result: 
Property: Let M be a real symmetric matrix defined as follows 

M=| m m„ |=(* M. (8) 
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Then we have 



V x |m\7 x (u(x,y)e z 



= -V ■ ( M _i (l(>t(M)V(/(.c. //) )e. 
Indeed, we note that 



V x |m\7 x (u(x,y)e z jj 

d / du du\ 
dx V dx dy) 

d ( du du\\ 

m iia e * 
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Furthermore, let M' be defined as 



m\ 1 m'i 



M' = ' T , 12 • (9) 

V. m 21 m 22 J 

We have 

V x (MV x r«(a;, y)e)j J = - V • I MVu J e 2 



J 2 ) 



if and only if 

. du . du du du 

m ll o 1" m l27T~ — m 22T: m ~^~ 

dx dy dx dy 



. du . du du du 

m 2iT- + m 22 — = mn- m— , 

dx dy dy dx 

which is true if M' = M -1 det(M). 

Using the above property, from ([7]), we derive that the transformed equation 
associated with (j3J) reads 

V • g^VH, + u 2 e^T~ l H z = , (10) 

with 

C} = e~ lf t /det(f ). (11) 



Here, T denotes the upper diagonal part of the transformation matrix T and 
T zz its third diagonal entry. 

Invoking the one-to-one correspondence (p2), we infer that the transformed 
equation associated with (jSj) reads 

V.^ 1 Vp + w 2 A- 1 r- 1 p = l (12) 

with 

£ T 1 = p 1 f/det(f). (13) 

In the sequel we will also consider a compound transformation. Let us 
consider three coordinate systems {u, v, w}, {X, Y, Z}, and {x, y, z} (possibly 
on different regions of spaces). The two successive changes of coordinates are 
given by the Jacobian matrices J x x and Jx u so that 

Jxu = JxxJxu ■ (14) 

This rule naturally applies for an arbitrary number of coordinate systems. 



4 Flying carpets over a flat ground plane 

This section is dedicated to the study of carpets levitating above a ground 
plane, that conceal to certain extent any object placed anywhere underneath 
them from plane waves incident from above. The construction of the carpet 
is a generalization of those considered in |39j to carpets flying over ground 
planes or located on either sides of rectangular objects, as shown in Fig. [1] 
with certain altitude y = y$. In the context of pressure waves, this corre- 
sponds for instance to the physical situation of a carpet which flies in mid-air 
if y = is the altitude of the ground, or a carpet which floats in mid-waters 
if y = stands for the bottom of the sea. This formalism allows us to study 
carpets which are either flying/floating on their own, or which are touching 
a cylindrical object on the ground or in mid-air/water. 

4.1 The construction of carpets over horizontal planes 

When the carpet is made on the y-axis (we call it a carpet flying above the 
x-axis), we consider a transformation mapping the region enclosed between 
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Figure 1: (Color online) Construction of a carpet above the a;-axis. The 
transformation ffl5|) shrinks the region between the two curves y = yo (dotted 
red) and y = y 2 (x) (dashed blue) and two vertical segments x = x and x = 
X\ into the region between the curves y = y\(x) (solid black) and y = y2{x) 
(dashed blue) and the vertical segments x = x and x — x\ (carpet). The 
curvilinear metric inside the carpet is described by the transformation matrix 
T, see (fT9|) . corresponding to density and bulk modulus of the metafluid 
given by (J6]). The grey rectangle can be either filled with air/ambient fluid 
(flying/floating carpet) or be replaced by a rigid cylinder. 



two curves (x,y ) and (x,y 2 (x)) to the one comprised between (x,yi(x)) and 
(x,y 2 (x)) as in Fig. [TJ where (x,y ) is mapped on (x,yi(x)) and (x, y 2 (x)) is 
fixed point-wise, of the form 



y> = a(x)y + (3(x) with a = » and (3 = ^gy 2 (15) 
z' = z 

So the inverse of this transformation is given by 

x = x' 

z = z' 

The above transformation ffTB"]) has the following Jacobian matrix 

J~* = !, {X ' y ' Z \ = \9 1 | . (17) 
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in which we have set 

dy 1 / d(3 , a\ da \ 
dx' a 2 \ dx' dx' J 

(V2 ~y') (V2 -yo) dy! 



(ya - yi) 2 dx 

(yi -y')(yi-yo)dy 2 
(2/2 - yi) 2 dx ' 



;i8) 



Hence we get 

T- 1 = J^JLjJdet^) = I -g (l + g 2 )a I . (19) 




4.2 Analysis of the metamaterial properties 

Let us now look at an interesting feature of the invisibility carpet. The 
eigenvalues of T _1 are given by: 

z = l,2 : Xi = — ( 1 + a 2 + g 2 a 2 + (-l)V-4a 2 + (1 + a 2 + g 2 a 2 ) 2 

A 3 = (20) 

a 

We note that Aj, i = 1,2, and A 3 are strictly positive functions as ob- 
viously 1 + a 2 + g 2 a 2 > \J— Aa 2 + (1 + a 2 + g 2 a 2 ) 2 and also a > 0. This 
establishes that T _1 is not a singular matrix for a two-dimensional carpet 
even in the case of curved ground planes, which is one of the main advantages 
of carpets over cloaks [39J. The broadband nature of such carpets remains 
to be investigated when one tries to mimic their ideal material parameters 
with structured media. 

The carpets in Fig. [2] and [3] are made of a semi-circle y = 2/2 (#) as an 
outer curve and a semi-ellipse y = yi(x) as its inner curve, both centered at 
(a , b\) with 



yo = bo, 3/1(2:) = bi + (1 - — )Jrl - (x - a ) 2 , 
y 2 {x) = b 1 + ^/r 2 -(x-a ) 2 . (21) 
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Plugging the numerical values r := 0.2, k = 0.1, a = in, one gets 
(5x(506i(6 - h) + 1 - 25x 2 + 50y(6i - M) 

# — 3 5 

(l-25s 2 )I 

where bo is the y-coordinate of the ground (typically, the ground can be taken 
as y — 0, so that 6 = 0) and &i is the hight, measured on the y-axis from 
the origin, at which the carpet is flying. 

4.3 The construction of carpets over vertical planes 

We note that for a carpet flying above the y-axis, if we consider the trans- 
formation mapping the region enclosed between two curves (xo(y),y) and 
(x2{y),y) to the one comprised between (x\(y),y) and (%2(y),y) as in Fig. [U 
where again (x (y),y) is mapped on (xi(y),y) and (x 2 (y),y) is fixed point- 
wise, that is, if 




a(y)x + j3(y) with a 

y 



z 




X2-XQ 



then the Jacobian of the inverse transformation now reads 




(22) 



where h is defined as 



h 




(x 2 — x')(x2 — Xq) dx 



+ 



(x 2 - xi) 2 dy 

(Xl — X'){X\ — Xq) dX2 



(23) 



(x 2 - Xi) 2 dy ' 



and hence we get 




(24) 



The same analysis as in Section fl~2l literally applies here, as well. 
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4.4 Numerical results for a carpet over a plane 

We first look at the case of a pressure plane wave incident upon a flat rigid 
ground plane (Neumann boundary conditions) and a carpet above it. The 
inner boundary of the carpet is rigid (Neumann boundary conditions). We 
report these results in Figure [2] where we can see that the altitude y' = 0.9 
leads to less scattering than the other two flying carpets. Of course, the 
carpet attached to the ground plane leads to perfect invisibility. 

We then look at the case of the rigid ground plane with a rigid circular 
obstacle on top of it. Some Neumann boundary conditions are set on the 
ground plane, the inner boundary of the carpet and the rigid obstacle. Once 
again, we can see in Figure [3] that the altitude y' = 0.9 for the flying carpet 
is the optimal one. 



(a) -1 (b) 




Figure 2: (Color online) 2D plot of the real part of the total pressure field 
dte{p): Scattering by a pressure plane wave of wavelength 0.15 incident from 
the top on a flat ground plane with a carpet above it, here yi(x) = b\ + 
1/2^0.04 - x 2 and y 2 (x) = bi + vCM" ) carpet touching the ground, 

(b) carpet flying at altitude b\ = 0.7, (c) at b\ = 0.9 and (d) at b\ = 1. The 
optimal altitude (for a flying carpet) b\ = 0.9 is noted. 
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Figure 3: (Color online) 2D plot of the real part of the total pressure field 
dte{p): Scattering by a pressure plane wave of wavelength 0.15 incident from 
above on a circular object of radius r = 0.6 lying on a flat ground plane; (a) 
obstacle without carpet, (b) with a carpet flying at altitude b\ = 0.7, (c) at 
b\ = 0.9 and (d) at b\ = 1. The optimal altitude b\ = 0.9 is noted. 



4.5 Numerical results for a carpet over a square cylin- 
der over a plane 

In this paragraph, we would like to render a pipeline lying on the bottom 
of the sea or floating in mid-water undetectable for a boat with a sonar at 
rest just above it on the surface of the sea. However, instead of reducing 
its scattering cross-section like in acoustic cloaks, we rather mimic that of 
another obstacle, say a square rigid cylinder. 

In Fig. HI we plot the total field when the pressure wave of wavelength 0.15 
is incident from above on a rigid square cylinder of sidelength 0.4. In panel 
(a) the obstacle touches the ground (e.g. bottom of the sea) at y — —1.26 
and is centered about x = 0. In panel (b), it has three carpets shaped as 
tents attached to its sides. The construction of the carpets is as follows: The 
right-most tent is defined by xo = 0.2, x\{y) = —y — 0.66, £2(2/) = —y ~ 0.56 
(upper part) and xi(y) = y + 1.46, x 2 (y) = y + 1.56 (lower part). The left- 
most tent is defined by x = —0.2, xi(y) = —y — 1.46, x 2 {y) = —y — 1.56 
(lower part), and Xi(y) = y + 0.66, x 2 (y) = y + 0.56 (upper part). The 
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pressure and bulk modulus of these left-most and right-most tents is deduced 
from the expression of the transformation matrix for vertical walls, see (I24p . 
Finally, the uppermost tent is defined by yo = —0.86, 2/1(2) = —x — 0.66, 
2/2(2) = —a; — 0.56 (right part) and 2/1(2) = x — 0.66, 2/2(2) = x — 0.56 (to the 
left). Here, we use the expression of the transformation matrix for horizontal 
walls, see f lT^l) . 

We now look at the case of a floating rigid square cylinder of sidelength 
0.4 in mid-water (e.g. a pipeline). In panel (c), the obstacle is flying on 
its own. In panel (d), it has three tentlike carpets on its sides and yet still 
scatters the incoming wave as in (c). The construction of the carpets in panel 
(d) is as follows: The right-most tent is defined by x = 0.2, Xi(y) = — y + 0A, 
22(2/) = —y + 0.5 (upper part) and Xi(y) = y + 0.4, x 2 (y) = y + 0.5 (lower 
part). Similarly, the left-most tent is defined by xo = —0.2, x\(y) = — y — 0.4, 
22(2/) = —y — 0.5 (lower part) and xi(y) = y — 0.4, 0:2(2/) = V ~ 0-5 (upper 
part). Finally, the uppermost tent is defined by yo = 0.2, 2/1(2) = — x + 0.4, 
2/2(2) = —x + 0.5 (right) and 2/1(2) = x + 0.4, 2/2(2) = x + 0.5(left). 

We note that while forward and backward scattering in panels (b) and 
(d) are not negligible (as would be the case for acoustic cloaks), these are 
instantly recognizable as that of a square obstacle on the ground (panel a) 
or floating in mid- water (panel c). The next question to address is whether 
such a generalized cloaking works at other incidences. 

4.6 Numerical results for a carpet over a square cylin- 
der over a plane in grazing incidence 

We now look at the case of grazing incidence, keeping otherwise the same 
configuration as in Fig. HI We note that the total field for the square obstacle 
is exactly the same with and without the three carpets, which is a further 
evidence that carpets allow for multi-incidence cloaking. Of course, one can 
hide any object (e.g. a semi-disc) under the carpets in panel (b) and this 
will scatter as a square. However, in the case of a flying object, it is required 
that the lower boundary of the hidden object be flat (e.g. a semi-disc) in 
order to mimic the diffraction pattern associated with a square obstacle. 

Taking into account that it should be possible to design broadband car- 
pets as their material parameters are non-singular, our numerics suggest that 
carpets are thus a very interesting alternative to invisibility cloaks. They can 
either reduce the scattering cross section of a rigid object or mimic that of 
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Figure 4: (Color online) 2D plot of the real part of the total pressure field 
dte{p): Scattering by a pressure plane wave of wavelength 0.15 incident from 
above on a rigid square cylinder of sidelength d = 0.4 lying on a flat ground 
plane on its own (a); with three (tentlike) carpets on it sides (b); flying over 
the ground plane on its own (c); with three (tentlike) carpets on its sides (d). 

another rigid object. 

Finally, we should emphasize that we numerically checked that our pro- 
posal for generalized carpets works fine for a pressure plane wave incident 
on the square at any angles ranging from normal to grazing incidence. This 
should not come as a surprise as the same holds true in the original proposal 
by Li and Pendry [39|. 

5 Flying carpets surrounding a circular cylin- 
der over a plane 

We propose in this section to design a carpet which geometrically looks like 
Pendry's cylindrical cloak (9), but in fact, it scatters any incoming plane 
wave just like a rigid cylinder of smaller radius Rq. The radius of such an 
effective rigid cylinder can be of any size smaller than that of the coated 
region. Such a carpet acts in a way similar to the invisibility cloaks proposed 
by Greenleaf et al. and Kohn et al. which are based upon the blow up 
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(a) 1 (b) 




1 



Figure 5: (Color online) 2D plot of the real part of the total pressure field 
dte{p): Scattering by a pressure plane wave of wavelength 0.15 incident from 
the right (grazing incidence) on a rigid square cylinder of sidelength d — OA 
lying on a flat ground plane on its own (a); with three tentlike carpets on 
it sides (b); flying over the ground plane on its own (c); with three tentlike 
carpets on its sides (d). 

of a small ball (of radius 77 -C 1) rather than a point, thereby leading to 
approximate invisibility. However, in our case, the radius of this ball (a disc 
in 2D) is finite, and our claim is that the carpet mimics the electromagnetic 
response of the rigid circular cylinder, just like previous carpets did for planes. 
The circular carpet actually plays the opposite role to the super-scatterer 
of Nicorovici, McPhedran and Milton whereby a cylinder surrounded by a 
coating of negative refractive index material scatters as a cylinder of diameter 
larger than the coating itself [5 J . 

However, whether this coated region is completely empty (e.g. Fig 0(b) 
) or many objects are hidden inside it (see e.g. Fig. [3(d) where the small 
cylinder of radius Ro has actually been put back inside the carpet, along with 
other objects), such a carpet produces a mirage effect that tricks an external 
observer into believing that this whole region is just the small rigid cylinder. 
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5.1 The construction of circular carpets 



The carpet consists of a cylindrical region C(0, R\) of radius R\ to be coated 
and the coating itself which is the space between an inner cylinder of radius 
Ri and an outer one of radius R 2 > R\. As above, the material properties 
of this coating will be deduced by pullback, via a transformation that fixes 
angles just like in [S]. The transformation in [5] is indeed the particular case, 
of the one under consideration here, associated with the value Rq = of the 
radius of an imaginary small cylinder as explained below. 

The carpet constructed below is a generalization of the flying carpet stud- 
ied in the previous section. The geometric transformation we consider here, 
is a smooth diffeomorphism of the closure of the outside M 3 \C(0, Rq) of a 
solid cylinder C(0,Rq) of radius -R , where Rq < R\. It coincides with the 
identity map outside the solid cylinder C(0, R2), fixing its boundary point- 
wise, but now maps the region A(Rq, R2) between the two coaxial cylinders 
of respective radii Rq and R2 into the space A(R±, R2) between the cylinders 
of radii Ri and R2, as in Fig. El More precisely, in A(Rq, R2) this geometric 
transformation can be expressed as 

r' = R 1 + a (r - R ) with a = g=g 

6' = 6 "" " (25) 

j „ 



with inverse 



z = z 

The Jacobian 3 rr ' of the latter is 



r = R + ±(r'-R 1 ) 

6 = 6' 
1 



d{r '°' Z) diag(I,l,l). (26) 



d(r',6',z') & V 

Let us denote by J xr the Jacobian of the change (r, 6, z) 1— ► (x, y, z) from 
Cartesian to polar coordinates and J rx := J~^. The Jacobian J xx i of the above 
transformation in Cartesian coordinates (x',y',z') 1— > (x,y,z) is obtained by 
applying the chain rule (JUJ) to get J xx i = J xr J rr >J r r x i, so that the tensor 
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3 XX ,J XX , det(J xx >) reads 



I l+(m 2 -l)cos 2 (6>) (m 2 -l)sin(6))cos(6)) g ^ 

m m 

-l)sin(fl)cos(fl) m 2 +(l-m 2 )cos 2 (e) q 

m m 

o o m 



V 

1 Til 

R{9) diag(m,-,— ) R(- 
m or 



(27) 



/cos(0) -sin(0) 0\ 
where R(#) = I sin(#) cos(#) is the matrix of the rotation with angle 

V l) 

9 in the x?/-plane and m = ^ = 1 — (f^rf^-fr- 




Figure 6: (Color online) Construction of a circular carpet of inner radius 
Ri (solid dark) and outer radius R 2 (dashed blue) from cylinder S (0, R ) of 
smaller radius Ro (dotted red). The transformation ( 125]) shrinks the whole 
hollow cylindrical region A(Rq, R2) of inner radius Rq and outer radius R2 
into its subset ^4(i?i,i2 2 ), the cylinder 5(0, Rq) being stretched to <S(0, 
whereas 5(0, R2) is fixed point- wise. Such a carpet, a metafluid with bulk 
modulus and density given by (127]) and ([6]), scatters pressure waves as the 
rigid cylinder «S(0, Rq). 
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Figure 7: (Color online) 2D plot of the real part of the total pressure field 
dte{p): Scattering by a plane wave of wavelength 0.15 incident from above 
on a flat ground plane and (a) a circular object of radius Ro = 0.2 flying 
at altitude on its own; (b) an empty cylindrical region of radius R\ = 0.5, 
the "coated" region, surrounded by a carpet of inner radius R\ and outer 
radius R2 = 0.7. This hollow cylindrical carpet is designed to behave exactly 
like the object in (a) alone, irrespective of the form of any other additional 
object that may be enclosed inside; (c) same object as in (a) with now three 
additional small rigid cylinders touching it. The scattering is clearly different 
from that in (a); (d) now the three objects in as (c) have been hidden inside 
the carpet (b) and yet an outer observer will not be able to tell the scattering 
in (a) and in (d) apart. 

5.2 Analysis of the metamaterial properties 

The tensor T -1 diagonalizes as diagfm, — , 2*) where m = 1 — .p 1 "^, ^j- is 
bounded from below and above as it satisfies < a^r < m < a. This means 

ill 

that the material parameters are not singular, unlike the cloaking case as in 

Obviously, one realizes that, when Rq tends to zero, one recovers the case 
in [9] where the material properties are no longer bounded, but one of them 
tends to zero whilst another one recedes to infinity, as we approach the inner 
boundary of the coated region, see also [13] . 
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5.3 Mirage effect for a cylinder surrounded by a carpet 



We report the results of our simulations in Fig. [7J for a circular pipeline 
which is floating in mid- water, see panel (a). We then replace this pipeline 
by a circular carpet, see panel (b), which reflects a pressure plane wave from 
above at wavelength 0.5 in exactly the same way. When we add three small 
pipelines to the original one, see panel (c), the reflected field is obviously 
much different. However, when we surround the four pipelines by the circular 
carpet, see panel (d), the reflected wave is that of the original pipeline. Such a 
mirage effect, whereby a rigid obstacle hides other ones in its neighborhood, 
can thus be used for 'sonar illusions'. For instance, an oil pipeline might 
reflect pressure waves like a coral barrier so that a sonar boat won't catch its 
presence. Unlike for earlier proposals of approximate cloaks [T3|[Ti] scattering 
waves like a small highly conducting object, we emphasize here that we start 
the construction of the circular carpet by a finite size disc. 

6 Multilayered circular carpet for broadband 
mirage effect 

6.1 Reduced material parameters for circular carpets 

We now want to simplify the expression of the inverse of the transformation 
matrix T in order to avoid a varying (scalar) density (resp. permeability 
in optics). For this, we introduce the reduced matrix T~ d = diag(a, -) 
which amounts to multiplying T -1 in (127)1 by a/m. We deduce from ([2]) the 
transformed governing equation for pressure waves in reduced coordinates: 



V r , e ■ diag(— , — )V r ,eP + ^ 2 A 1 ap = , 



a a 



(28) 



and similarly for transverse electric waves using ([3]): 



V r . 9 ■ diag( 



a 



a 



)V r ,eH z + u 2 e fi aH z = , 



(29) 




r 2 -Ri 
R2-R0 ' 
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6.2 Homogenized governing equations for optics and 
acoustics 

To illustrate our paper with a practical example, we finally choose to design 
a circular carpet using 40 layers of isotropic homogeneous fluids. These 
fluids have constant bulk modulus and varying density. We report these 
computations in figure [HJ 

It is indeed well known that the homogenized acoustic equation for such 
a configuration takes the following form: 

V r , .(po 1 g l ^r,eH z ) + u 2 < X- 1 > p = , (30) 

where < A -1 >= J Q A -1 (r) dr and with p a homogenized rank-2 diagonal 
tensor (an anisotropic density) p = Diag(p r , pg) given by 

£ = Diag(<p- 1 >~\<p>) . (31) 

We note that if the cloak consists of an alternation of two homogeneous 
isotropic layers of fluids of thicknesses g?a and ds, with bulk moduli A^ and 
Xb and densities pa and ps, we have 

1 1 ( l_ + rj_ 



p r 1 + 7] \p A p B 

_ PA + VPB >_i 1 / 1 V 

Pe - — — , < A >- — — h — 

1 + T] 1 + 7] \X A X B 

where rj = ds/dA is the ratio of thicknesses for layers A and B and dA + ds = 
1. 

Using the change of variables 

Vrfi.ie^g 1 V r , e p) + e p uj 2 H z = 0, (32) 

where e is a homogenized rank-2 diagonal tensor (an anisotropic permittivity) 
g = Diag(£ r , e B ) given by 

e = Diag(< e~ l > _1 , < e >) . (33) 
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6.3 Acoustic paradigm: Reduced scattering with larger 
scatterer 

The acoustic parameters of the proposed layered circular carpet are therefore 
characterized by a spatially varying scalar bulk modulus p and a spatially 
varying rank 2 density tensor p given by (|33|) . We can further simplify the 
problem by choosing reduced acoustic parameters, so that the bulk modulus 
A is now constant, and all the variation is reported on the density, see figure 
[HI More precisely, pa varies in the range [0.1890; 0.5493] and ps varies in 
the range [1.7987; 2.1472]. We checked that this carpet is broadband as it 
works over the range of wavelengths A G [0.2, 1.4286], see Fig. [8] and Fig. [9j 
a multilayered carpet of radius 1 surrounding a rigid obstacle of radius 0.32 
scatters waves just like a rigid obstacle of radius Rq = 0.2. We note that the 
lower bound for the range of working wavelengths corresponds to the rigid 
obstacle we want to mimic. 

7 Conclusion 

In this paper, we have proposed some models of flying carpets which levitate 
(or float) in mid-air (or mid- water). Such cloaks can be built from acoustic 
metafluids: as explained by Pendry and Li in a recent work, one can for 
instance emulate required anisotropic density and heterogeneous bulk mod- 
ulus with arrays of rigid plates with a hemispherical sack of gaz attached 
to them [12] • But other designs proposed by Torrent and Sanchez-Dehesa 
would work equally well [17]. However, such flying carpets lead to some ap- 
proximate cloaking as they do not touch the ground (the inner boundary of 
the carpet in the original design of Pendry and Li is attached to the ground 
[5§] . Interestingly, other authors also looked at quasi-cloaking with simplified 
carpets [43] . We have also explained how one can hide an object located in 
the close neighborhood of a rigid circular cylinder, which in some sense can 
be classified as an external cloaking whereby a large scatterer hides smaller 
ones located nearby. Such an ostrich effect (which buries its head in the 
sand) has already been observed in the context of dipoles and even finite size 
obstacles located closeby a cylindrical perfect lens which displays anomalous 
resonances [6], [H] . However, here the coating does not contain any negatively 
refracting material, and this is an anisotropy-led rather than plasmonic-type 
cloaking mechanism. Actually, it is possible to use complementary media to 
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Figure 8: (Color online) 2D plot of the real part of the total pressure field 
dte{p): Scattering by a plane wave of wavelength 0.2 incident from above 
on a flat ground plane and (a) a circular object of radius Ro = 0.2 flying 
at altitude on its own; (b) an empty cylindrical region of radius R\ = 0.32, 
the "coated" region consisting of 40 layers of isotropic homogeneous fluid, 
see closer view in (c), of constant bulk modulus and density p, given in (d), 
surrounded by a carpet of inner radius R\ and outer radius R 2 — 1. The red 
curves represents the variation of pe = m 2 /a with respect to r G [0.32; 1]. 
The piecewise constant blue curve is a staircase approximation of the red 
curve, considering an alternation of 40 layers of density pa € [0.1890; 0.5493] 
and pb G [1.7987; 2.1472] using the homogenized formula (155)1 . 

cloak finite size objects (rather than only dipoles) at a finite distance |45) . 

We have discussed some applications, with the sonar boats or radars 
cases as typical examples. Another possible application would be protecting 
parabolic antennas from the negative impact of their 'supporting cable'. The 
feasibility of such carpets is demonstrated using a homogenization approach 
enabling us to design a multi-layered acoustic metafluid leading to a mirage 
effect over a finite range of wavelengths. We note that the route towards 
anamorphism discussed in this paper is very different from the proposal of 
Nicolet et al. [46J whereby an anisotropic heterogeneous object is placed 
within the coating of a singular cloak to mimic the scattering of another 
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(a) (b) 




(c) (d) 




Figure 9: (Color online) 2D plot of the real part of the total pressure field 
3?e(p) for the same geometric and acoustic parameters as in Fig. [8j Scattering 
by a plane wave of wavelength 0.5 (a-b) and 1.4286 (c-d), incident from above. 

object. 

Moreover, all computations hold for electromagnetic carpets built with 
heterogeneous anisotropic permittivity and scalar permittivity, which could 
be emulated using tapered waveguides as in [17]. We therefore believe that 
the designs we proposed in this paper might foster experimental efforts in 
approximate cloaking for both acoustic and electromagnetic waves. 
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